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A Sum-Product Model as a 
Physical Basis for Shadow Fading 
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Shadow fading (slow fading) effects play a central role in mobile communication system design 
and analysis. Experimental evidence indicates that shadow fading exhibits log-normal power distribution 
' almost universally, and yet it is still not well understood what causes this. In this paper, we propose a 

versatile sum-product signal model as a physical basis for shadow fading. Simulation results imply that 

o 

the proposed model results in log-normally distributed local mean power regardless of the distributions 
of the interactions in the radio channel, and hence it is capable of explaining the log-normality in a 
wide variety of propagation scenarios. The sum-product model also includes as its special cases the 

ON ■ 

conventional product model as well as the recently proposed sum model, and improves upon these by: 

(N ; 

a) being applicable in both global and local distance scales; b) being more plausible from physical point 
of view; c) providing better goodness-of-fit to log-normal distribution than either of these models. 

O 

Index Terms 



X 

Lognormal, shadow fading, channel modelling. 



I. Introduction 

In physical sciences, a theory of an observable phenomenon is developed by comparing 
theoretical predictions to measurements, and, possibly, modifying the theory if disagreement 
is found. In wireless radio channel research, one attempts to deterministically or stochastically 
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predict the behavior of radio waves in real-world environments by comparing channel models 
to measurements. If the model disagrees with experimental data, it should be revised. 

The purpose of this paper is to propose an improved model for the so-called shadow fading 
phenomenon observed in mobile radio channel measurements. Shadow fading (slow fading, 
shadowing) is loosely defined as the fluctuation in the received power averaged over a small 
area, typically having a diameter of 10 — 40 wavelength in outdoor environment [1], [2]. It plays 
a central role in various communication system design tasks, including coverage prediction for 
network planning and analysis of system-level radio network algorithms. In a large number of 
measurements conducted worldwide, the local mean power has been found to be well approx- 
imated by the log-normal distribution over a wide range of propagation environments, carrier 
frequencies, and transmitter-receiver distances [l]-[5]. Yet, no generally accepted explanation 
for shadow fading or its log-normality has been presented. Considering the great variety of 
propagation scenarios exhibiting log-normal power variation it seems clear that the central limit 
theorem, in one form or another, masks away the intricacies of the highly complex underlying 
propagation physics. Otherwise, it is difficult to explain why the log-normal statistics occur in 
environments where the dominating propagation mechanisms are different; for example, in urban 
and suburban macro cells [6], [7] and urban micro cells [8]. 

The conventional product model assumes that the radio signal propagating from the transmitter 
to the receiver experiences a cascade of random attenuations. From the central limit theorem 
it follows that the logarithm of the random product is approximately normally distributed [1], 
[3]. While this model appears physically plausible, experimental observations disagree with it to 
some extent. First, results from extensive computer simulations [9], [10] imply that the number 

og-normal distribution to 
is possible, the resulting 
path loss would be very high. Furthermore, measurements in urban areas indicate that most of 
the signal power is transported by propagation modes that experience much less than 10 channel 
interactions [11]. Second, the product model implies that path loss increases exponentially with 

'in this paper, by 'interaction' we mean either reflection, diffraction, or scattering. 



of random attenuations needs to be quite large for the convergence to 
take place, e.g. 30 — 100. While such a large number of interactions^ 
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distance d, and that the standard deviation of shadow fading is also distance-dependent [12]— [15] . 
However, majority of measurements suggest a oc d a path loss law (typically with a = 3 ... 7), 
and that the standard deviation of shadow fading is approximately independent of the receiver- 
transmitter separation [4], [5]. Hence, even if the product model is true, it seems that the number 
of measurable interactions is limited, so that e.g. exponential path loss difficult to detect in 
practiceo Third, an interesting measurement result with two receive antennas [18] indicates that 
when the patterns of the antennas are not identical (i.e., the antennas see different multipath 
environment), the correlation coefficient between the shadow fading processes observed at the 
antenna outputs was clearly less than unity. The result contradicts the product model which 
implicitly assumes that all multipaths experience the same multiplicative attenuative process. 
The observation that the variation of the local mean power is at least partly related to multipath 
propagation (instead of 'shadowing' by obstacles) is also supported by [19]— [21]. 

The sum model proposed in [14], [19] is based on the idea that the amplitudes of the planewaves 
impinging at the moving receiver change slowly in time/space with respect to small-scale fading, 
giving rise to slow variation in the local mean power. In a line-of- sight microcell this was shown 
to result in Nakagami-m or log-normal distributed slow fading [19]. The idea was generalized to 
general environments in [14], where the variation of the power of the plane waves was modelled 
stochastically, and it was argued based on the central limit theorem and a simple linearization 
approach that the resulting local mean power may be perceived as log-normal, regardless of the 
power distributions of the individual plane waves. Unfortunately, the applicability of the sum 
model in [14] is limited to a small area, called "extended local area", and it cannot explain the 
measured 5 — 12 dB path loss standard deviations in a larger, global area. 

It is evident that neither the product model nor the sum model can alone explain all aspects 
of the shadow fading phenomenon, and its log-normality in particular. The contribution of this 
paper is to introduce a unified sum-product model that has the advantage of being more plausible 
from physical propagation point of view than the product model, as well as being applicable in 

2 Some evidence to indicate that path loss is exponential has been reported in [16], while measurement results in [17] suggest 
that standard deviation of shadow fading is distance dependent. 
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a larger geographical area than the sum model. It is shown that, in addition to providing a more 
realistic physical basis, the shadow fading distribution induced by the sum-product model gives 
better fit to the log-normal distribution than the product model. Moreover, unlike with the sum 
model [14], the standard deviation of shadow fading can assume large values often encountered 
in measurements (e.g. 8—12 dB), since the justification of log-normality is not based on a linear 
approximation of the logarithm function, which only holds if the deviation of shadow fading 
relative to its mean is small. 

As in [14], we adopt a stochastic channel modelling approach in the sense that we do 
not specify the underlying physical propagation mechanisms (diffraction, reflection, scattering). 
Instead, the channel interactions are modelled probabilistically. This approach has the advantage 
that the results hold more generally, rather than for any given propagation environment. 

The paper is organized as follows. In Section [TT] we introduce the sum-product model, and in 
Section [HI] it is shown how the product and the sum models are special cases of it. Simulation 
set-up is detailed in Section [IV] while the simulation results confirming the applicability of the 
sum-product model are given [V] Section [VI] concludes the paper. 

Notation: A scalar, vector, and matrix are denoted by a, a, and A, respectively, while A T 
and A^ stand for the transpose and conjugate transpose of A. The expected value of a random 
variable X is denoted by E[X], A set {xi,x 2 , . . . ,xn} is also denoted with {x n }^ =1 , or just 
{x n } when N is obvious from the context. 

II. Sum-product model for radio propagation 

As the first step of this section we explain the fully determistic point-to-point version of the 
sum-product model, where all channel interactions are assumed constant (transmitter, receiver 
and channel are not moving). In the second phase, we let the receiver move within a small local 
area, such that the phases of the received plane waves vary in a random fashion, and we then 
examine the signal power averaged over small-scale fading. In the third step, the receiver is 
allowed to move in a large, global area, and all channel interactions are assumed random. 
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A. Point-to-point response 

We introduce the sum-product model by means of a simple example. Consider the point-to- 
point radio propagation scenario shown in Fig. [TJ Crosses in the figure represent interactions 
experienced by the propagating plane waves (or, "rays"). The interactions in the radio channel 
can be either scattering, reflection, or diffraction. The rays with complex responses {bi, b 2 } first 
interact with the first layer of interacting objects, and then continue onwards to interact with the 
second layer, and so on. At the receiver antenna output, the propagating rays are weighted by 
{a\, a 2 } and summed to yield a complex scalar received signal y. Some remarks on the model: 

• It is assumed that {a n }, {b m }, and {sk,ij} include the directional responses of the interacting 
objects at both ends of the connecting edge, including the free space path loss in between 
them. 

• It is assumed that all channel interactions have directional responses. To give an example, 
consider a channel interaction (e.g. wall reflection) excited by a plane wave from a direction 
6 1. We then assume that the response of the channel interaction towards another direction 
6*2 can be written as a product of two directional responses Sfc_i(0i)sfc(0 2 )- 

• Example 1 (see Fig. [7]): the complex number b\ captures the complex responses of the 
transmitter antenna to the direction of the channel interaction l A\ radio channel between 
antenna and l A\ and the channel interaction 'A' towards the direction of the transmitter. 

• Example 2 (see Fig. [7]): the complex number {si j2 i} captures the complex response of the 
channel interaction 'A' towards the direction of channel interaction channel response 
between interactions 'A' and '£>', and response of channel interaction '£>' towards channel 
interaction 'A'. 

This received signal at a fixed transmitter and receiver position can be written as 



V = dl [S 2 ,ll0l,llfrl + ^1,1262) + S2,12(S1,21&1 + Sl,22&2)] 

+ a 2 [s2,2l(si,ll&l + Si, 1262) + S 2 ,22(S1,21&1 + S 1}Z iki)] , 

or in matrix form 
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y = a T S 2 Sxb, (1) 

where a = [oi a 2 ] T , b = [bi b 2 ] T , and 



k = {l,2}. (2) 



Sfc,21 s k,22 

All scalars {a n }, {6 m }, {s&,y } are assumed complex valued. For now, we assume the propagation 
environment can be modelled as a network of channel interactions, as shown in Fig. [H where 
direct connectivity exists only between interactions of neighboring "layers". For example, the 
interactions at the transmitter end in Fig. Q] are connected to the interactions at the receiver end 
via the two intermediate ones. Hence the graph is not fully connected, which, interpreted in 
terms of radio propagation, means that the path loss between certain channel interactions is so 
large that it can neglected. In the sequel, we shall call a vertical group of interactions a "layer 
of interactions". Hence, in Fig. [H there are three layers of interactions, with two interactions 
in each. A matrix characterizing the interactions between two neighboring layers is called a 
"coupling matrix". 

More generally, we assume the signal model depicted in Fig. |2l which suggests that the 
received signal can be written as a bilinear form 

y = a^Sb, (3) 

c 

= a T c, (4) 

where 



a = 


[at ■■■ 


a N f 


b = 


[h ••■ 


b M ] T 


c = 


[ci ••• 


c N f 


S = 


SkSk- 


-l ■ • • 



(5) 
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' a y er ° f interaction 




Fig. 1. Example of the sum-product signal model with K — M — N — 2 (three layers of channel interactions). 

The composite coupling matrix S characterizes how the transmitted plane waves are coupled to 
the received ones. In general, {Sfc}^L 1 need not be square. By allowing rectangular matrices, one 
can vary the number of interactions within a layer. Empty interaction corresponds to a column 
of zeros in Of course, for © to make sense it is required that Sk has N rows and that 
Si has M columns, so that S is N x M. By imposing structure on the coupling matrices the 
bilinear form © lends itself for modelling a great variety of radio propagation scenarios; some 
examples will be given in Section [Till 

B. Signal power averaged over small-scale fading 

While the model in © may be applicable to wider range of purposes, in this paper we use 
it to study the distribution of received signal power averaged over small-scale fading. The local 
mean power is obtained by averaging over small-scale fading, or, equivalently, averaging over 
the phase of the received signal y. Towards this end, suppose now that the receiver in Fig. |2] is 
moving. To average over the small-scale fading we assume that only the phases of the entries of 



February 1, 2008 



DRAFT 



s 




transmitter 



receiver 



equivalent 
coupling matrix 



Fig. 2. Signal model. S is the N x M composite coupling matrix describing the mapping between the transmitted and received 
plane waves. 



a are changing independently and uniformly over [0, 2ir). All other quantities, including {|a n |}, 
{b m } and {s^}, are assumed constant within the area over which the average is taken[]. In 
other words, the small-scale fluctuation of y is assumed to be entirely due to variation in the 
phases of {a n }. If the aforementioned assumptions hold, the local mean power is obtained by 
taking the expected value over the phases of a: 

P = EN 2 ] 

= E[|a T Sb| 2 ] 

= E[b^S H aa H Sb] 

= b H S"E[aa H ]Sb 

= b H s H r a sb, (6) 

where S was defined in ©, and T a = E[aa H ] is an iV x N diagonal matrix with {|a n | 2 }^ =1 on 
the diagonal. Alternatively, we may write © as 

3 In practice, one would approximate the theoretical mean power with e.g. sliding mean over measured power of the small-scale 
fading signal. In measurement analysis, the averaging length is conventionally chosen as about 10 — 40 wavelengths in outdoor 
radio propagation environments [1]. 
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N 

P = Y,W?\Cn\\ (7) 
n=l 

where {c n } denote the entries of c [see ©]. This expression has the same form as the one 
given for local area power in [14], although it has been derived from the viewpoint of the more 
general signal model considered in this paper. 

C. Variability of local area power (P) - shadow fading 

In CD), it is assumed that {|a n |} and {c n } are constant, and hence P is constant for a given 
averaging area, also called local area [14], [22]. Consider now an experiment where the power 
is measured over a large number of local areas, as shown in Fig. |3] Assuming that the number 
of rays at the transmitter (M) and receiver (iV) end as well as the number of "layers" (K) in the 
signal model © are constant over all measured local areas, i.e., for all indices q in Fig. [3l we 
obtain a set of local area powers: {P q }® =1 . By further assuming that {a n }, {b m }, {sk,ij} for all 
n, m, k, i,j can be modelled as random variables, we can treat {P q } as realizations of a random 
variable with an unknown distribution. From the above assumptions, it follows that the average 
local area power, which can be interpreted as average path loss, is constant. This can be seen 
by taking the expected value of ©, where {|a n |} and {c n } are by assumption stationary. 

In this paper, we are interested in the distribution of the local area power, P. In Section [V] 
we provide simulation results that imply that the sum-product model © results in log-normally 
distributed local area power, regardless of the distributions of the individual components of the 
model. Fig. |4]gives a preview of this interesting result. The empirical cdf of 10 log 10 (F) in Fig. @] 
is from a Monte Carlo simulation where {|a n |}^.=i> {|^m|}m=i' an( ^ {l s Mjl}fc=i are from uniform 
and R distribution, to be defined in Section [IV] Normal distribution is plotted for reference with 
dotted line. It can be seen that the empirical distribution of 101og 10 (P) is very close to normal, 
i.e., P is log-normally distributed. 

Note that by setting N = M = 1 and assuming that the coupling matrices are scalars, 
© collapses to a product of K + 2 scalar random variables, which is known to converge in 
distribution to log-normal under mild conditions as K — ► oo. Remarkably, the simulation results 
in Fig. |H and Section |V] suggest that the central limit theorem for the logarithm of a product of 
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local area 




Fig. 3. Sampling of mean power over several local areas each having a random set of parameters {a q }^ =1 , {c 9 }^ =1 . In Section 
III-BI only {a 9 }^ =1 are assumed random, while in Section Hl-CI also {c 9 }^ =1 are random. 



random positive scalars generalizes to the matrix quadratic form in ©; this is itself an interesting 
observation that requires further investigation. 

III. Special cases of the general signal model 

Before simulation results, we briefly illustrate how the general sum-product model yields as its 
special cases the classical multiplicative shadow fading model and the recently proposed additive 
model. 

A. Multiplicative global shadow fading model 

In the classical multiplicative shadow fading, all plane waves experience a cascade of random 
multiplicative attenuations. This can be modelled by settingj S k = s k l N in ©. In this special 



4 We assume for simplicity that all Sfc are square. 
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Probability plot for Normal distribution, K=40, N=M=10 




-50 -40 -30 -20 -10 10 20 30 40 

data [dB] 



Fig. 4. Example of empirical distribution of 101og 10 (P) with P given in $6$. The amplitudes of the complex entries of a, b, 
and {Sfej^lj are independently drawn from R and uniform distribution over the unit interval, and the phases are uniform over 
[0, 2ir). Normal distribution is shown (dotted line) for reference. The R distribution is defined in Table U Section ITVl 



case, © reduces to 



N \ K 



P= Z>n| 2 N 2 -\{\Sk\ 2 . (8) 

\n=l / k=l 

This is a product of K + 1 random factors, whose distribution is well-known to converge to 
log-normal for large K. 



B. Additive local shadow fading model 

In [14], an additive model was proposed as a physical basis for shadow fading with an area that 
consists of a several local areas. From the viewpoint of the present framework, this "extended 
local area" is on one hand assumed large enough so that {|a n | 2 } are randomly time-variant, 
and on the other hand small enough so that c, i.e., the intermediate propagation process, can 
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be assumed constant. If these conditions hold, {|c„| 2 } are constant in ©, and shadow fading is 
only due to variation of the amplitudes of the plane waves at the receiver end. It was shown in 
[14] that the resulting distribution of P may be perceived as log-normal if the standard deviation 
of P relative to its mean is not too large, or if {|a n | 2 }^ =1 themselves are log-normal and N is 
small. 

C. Other special cases 

By suitable selection of the coupling matrices {S^}, it is possible to model a variety of 
propagation scenarios. We provide three examples. 

1) Line-of- sight: In line-of-sight propagation, the direct ray experiences only free space loss, 
while other rays interact with the environment. This scenario can be modelled, for example, by 
setting 



In the above pi denotes the free space path loss, nX m an n x m matrix of zeros, and Sfc in i os 
the kth coupling matrix of the non-line-of- sight rays. We have assumed for simplicity that all 
coupling matrices are square. 

2) Keyhole propagation: It has been shown theoretically that a keyhole degrades the infor- 
mation capacity of a MIMO communication system [23]. Although in this paper we focus on 
the single-input single-output case, we remark that the signal model © could also be extended 
to MIMO case. To model a keyhole, one would simply set one of the coupling matrices to a 
row vector. 

3) Cluster-based propagation: Consider the case where the plane waves propagate in L 
clusters so that each cluster of rays is completely decoupled from the other clusters. This scenario 
is realized by letting {S^} be block diagonal, i.e., 




S*,i o 



(9) 







s 



k,L 
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where {S^ denotes the A;th coupling matrix of the /th cluster. Further simplification results 
by assuming that all rays within a cluster undergo the same multiplicative attenuation, which 
amounts to letting = s^I, where I is an identity matrix. 

IV. Simulation procedure 

In order to examine the distribution of shadow fading, we can let b, {S^}, and a (also called 
"components" in the sequel) in the quadratic form © be random variables. Unfortunately, we 
have been unable to find in statistics literature any analytical results or limit theorems on the 
distribution of the quadratic form. Therefore, in the sequel we shall resort to Monte Carlo 
simulations in order to study its distribution empirically. The empirical findings, reported in 
Section |Vj indicate that the distribution of the quadratic form converges to log-normal, regardless 
of the distributions of its components. 

A. General assumptions 

The following assumptions are made to restrict the number of free parameters in the simu- 
lations. Computer experiments indicate, however, that the conclusions in Section [V] hold also 
under more general conditions. 

• We set M = N and assume that the coupling matrices {Sk}^ =1 are all square matrices. 

• It is assumed that phases of {a n }, {b n }, and are independent and uniformly dis- 
tributed over [0,27r) for all n,k,i,j. It is also assumed that the phases are distributed 
independently of the corresponding amplitudes {|a n |}, an ^ {| s fc,i?'l}- 

• It is assumed that the amplitudes {|a n |}, {|6 n |}> and {\s k are identically and indepen- 
dently distributed for all n,k,i,j. 

In summary, the tunable parameters in the Monte Carlo simulation are: the number of plane 
waves (N), the number of coupling matrices (K), and the amplitude distribution of the compo- 
nents. 

B. Distribution of the components 

To our knowledge, there are no published measurement studies available on the distributions 
of the attenuations of the channel interactions. To keep the amplitude distributions physically 
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reasonable, the support of the random variable should be restricted to unit interval [0, 1], since 
the radio channel is a passive component and hence the energy of the propagating signal 
cannot increase in interactions. Assuming a positive random variable X defined over [0, oo], 
the transformation Y = (1 + X)" 1 , for example, results in another random variable Y defined 
over the unit interval. For purposes of simulations, we shall assume that X is either Rayleigh, 
with density 



2 

X X 



fx(x) = — exp f -— ) , (10) 



or log-normal with density 



The corresponding distributions of Y will be called R and L distributions, for short. In addition, 
we shall employ the beta distribution, which for A = B = 1 reduces to the uniform distribu- 
tion. The amplitude distributions are summarized in Table HI where the same notation as with 
MATLAB's random number generators has been adopted. The distributions of the attenuations 
(in dB scale) are illustrated for selected values of parameters in Fig. |5J 

To further keep the simulation model physical, one should normalize the columns of each 
coupling matrix by \/N, or 



where the entries of Sk are from a distribution defined over the unit interval. This normalization 
ensures that the propagating energy is not amplified in channel interactions. However, it is easy 
to see by plugging (fl2l) into © that this normalization only shifts the mean of the distribution of 
In P but does not affect its shape. As we are primarily interested in the shape of the distribution 
of the random variable In P, we choose to ignore the power normalization (fl2l) in the simulations. 
Furthermore, for easier visual comparison all distributions shown will also centered so that the 
mean value is dB. 
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TABLE I 

Amplitude distributions used in the Monte Carlo simulations 



name 


definition^ 


parameters 


beta 


Pdf: ^V^-V) 3 - 1 


{A,B} 


R 


Y — , X ~ Rayleigh 


B 


L 


Y = jtx > ~ log-normal 





' 13(^4, £?) denotes the beta function. 



C. Quantifying goodness-of-fit 

As in [9], to quantify the distance of the empirical cdf of In P from the normal distribution, we 
utilize the Kolmogorov-Smirnov (K-S) test statistic [24]. The test statistic is defined as the largest 
absolute difference between the empirical cdf (from simulation) and the normal distribution. Note 
that, again as in [9], our main interest is not to test the hypothesis that the observed data comes 
from the log-normal distribution^, but rather to compare the relative rates of convergence with 
different parameters and component distributions. This will allow us to make some conclusions 
of how well the sum-product model can explain the shadow fading phenomenon. 

V. Simulation results 
The purpose of the numerical experiments in this section is to answer the following questions: 

• How does the relative convergence speed of the sum-product model and the conventional 
product model depend on the number of interaction layers (K)1 What are the standard 
deviations of the resulting distributions, and how sensitive they are to the value of Kl 

• What is the impact of the number of plane waves (AT)? 

• What is the impact of the distribution of the attenuations of the channel interactions? 

In all cases, the empirical cdf and the corresponding K-S test statistic has been computed 
from a Monte Carlo simulation with 10 5 realizations of the random variable P [cf. ©]. 

5 For any nontrivial values of the parameters N,M,K and any distributions of the entries of {Sfc},a, b in ® it ls highly 
unlikely that the true distribution (not known in analytic form) is exactly log-normal. Therefore, testing hypothesis that InP 
is log-normal does not make sense since for large sample sizes the test will anyway fail with high probability. We conjecture, 
however, that as K —* oo, In P (suitably normalized) converges in distribution to normal. 
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-50 -40 -30 -20 -10 

t[dB] 

Fig. 5. Example of cdfs. 

A. Sum-product model versus product model 

The main focus of this paper is to study whether shadow fading could result from the sum- 
product signal model, and whether the results predicted by it are more plausible than those 
obtained from the 'text book' product model. To shed some light on these issues, we plot the 
value of K-S test statistics for increasing number of multiplicative layers, K. The test is repeated 
both for the product model in ® and the sum-product model in ©. Fig. [6] shows the results 
for the case where the amplitudes of the entries of {S^}, a, and c are independently generated 
from beta, R and L distributions (see Table HI) with parameters A = B = 1 (uniform), B = 10, 
and [i — a — 1, respectively. The number of rays is N = M = 10. 

From Fig. |6] it can be seen that for a given value of K, the empirical power distribution 
generated by the sum-product model is much more similar to the log-normal than that of the 
conventional product model. As implied by extensive numerical studies [9], [10], the convergence 
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of the product model requires a large number of factors (e.g. > 30) to converge to log-normal 
distribution. On the other hand, the convergence of the sum-product model is quite fast for all 
tested distributions, and the fit is also considerably better. 

When comparing models one should, in addition to the shape of the signal power distribution, 
also consider the mean path loss and the magnitude of its standard deviation. Unfortunately, it 
is not straightforward to compare models in terms of absolute path loss without incorporating 
some propagation physics into the model. This, however, would make simulations hopelessly 
complicated^], and therefore in this paper we limit our focus to the stochastic approach. Nonethe- 
less, it is of interest to compare the standard deviations of the distributions in Fig. [6] for different 
values of K, see Table lU From the table, it can be seen that the standard deviation with the 
product model increases quickly with K. With the uniformly distributed channel interactions, 
the deviation is 55 dB for K = 40, whereas for the sum-product model it is only about 9 dB. 
In general, with the sum-product model the standard deviation is less sensitive to K; for K 
increasing from 5 to 20 the value of the deviation ranges from 3.8 to 8.9 dB; these values are 
well within the range of typical values encountered in measurements. From radio propagation 
viewpoint, the value of K cannot be large, since this would result in very low received signal 
power owing to the multiple layers of interactions. In this sense, the relative insensitivity of 
the sum-product model to the value K can be related to the empirical observation that shadow 
fading standard deviation is relatively independent of the receiver-transmitter distance [3]. On 
the other hand, the standard deviation of the product model is quite sensitive to the number of 
attenuations; little or no empirical evidence of such behavior has been observed in measurements. 

B. Impact of number of planewaves 

Fig. [8] illustrates the impact of the number of planewaves (N) to the distribution of the sum- 
product model. Changing iV has only a minor effect on the distribution of the product model. 
This is because in © N affects only one term in the product, whereas in the sum-product 
model changing iV changes also the dimensions of the coupling matrices (since we assume for 

6 Numerical experiments with the simple product model incorporating over-the-rooftop propagation effects was considered in 
[10]. Similar experiments with the sum-product model would be even more complicated. 
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Fig. 6. Sum-product model versus product model. Small value of the test statistics indicates good fit. Critical value of the 
Kolmogorov-Smirnov test for 95% significance and 10 3 samples is shown for reference. Beta: A=B=1. R: B=10. L: /i = a = 1. 



TABLE II 
Standard deviations for Fig. [6] 



model/pdf 


std for K = 1/5/10/20/40 [dB] 


sumprod/beta 


2.7/3.8/4.9/6.6/9.1 


sumprod/R 


4.2/5.6/6.9/8.9/12.0 


sumprod/L 


3.1/4.2/5.3/7.0/9.5 


prod/beta 


9.0/19.5/27.5/38.8/55.1 


prod/R 


6.1/11.4/15.6/21.7/30.6 


prod/L 


6.7/14.1/19.6/27.7/39.0 
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Normal Probability Plot, K=5, N=M=10 



0.9999 








0.999 








0.99 






/ / ■ / J ' 


0.95 








0.9 








0.75 








0.5 








0.25 




product (beta) 




0.1 








0.05 




product (R) 




0.01 








0.005 








0.001 
0.0005 




i .'i 


: : \ sum-product (beta) 

;) 1 : sum-product (R) 

.' i ■ i i 



-80 -60 -40 -20 20 

data [dB] 



Fig. 7. Examples of empirical distributions of 10 log 10 (P) from Fig. [6]for K = 5. The straight dotted line displays the normal 
distribution. See Fig. g] for K = 40. 



simplicity that they are N x N). With all tested distributions, the sum-product model gives a 
better fit for all values of N. Interestingly, with the sum-product model the fit improves as N 
increases for the uniform and L distributions (for the R distribution no conclusions can be made 
about convergence for the examined range of N). Such convergence to log-normal distribution 
is somewhat surprising, since one would expect exactly the opposite since increasing the length 
of a and b implies increased averaging (summation) in the signal model. It appears that the 
increase in the size of the N x N coupling matrices balances the averaging effect. 

The averaging effect from increasing iV can be seen in the standard deviation of the sum- 
product distribution, see Table Hill The standard deviation of the distribution decreases sharply 
as iV increases. With the product model, TV has little effect on the standard deviation, which is 
in line with the earlier observation that number of rays affects only the distribution of the first 
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Fig. 8. Impact of number of rays, N, for different distributions. Small value of the test statistics indicates good fit. Critical 
value of the Kolmogorov-Smirnov test for 95% significance and 10 3 samples is shown for reference. 



factor in ©; this will not have much effect unless K is very small. 

C. Impact of distributions 

Fig. [9] illustrates the impact of distribution parameters on converge speed. Some remarks on 
the results: 

• With all tested interaction distributions, including many others not shown in Fig. |9l exhibit 
convergence to the normal distribution. This implies that (j6]) converges to log-normal under 
very mild statistical assumptions on elements of the component matrices and vectors. As we 
are not aware of any limit theorems that would apply to the signal model ©, the requied 
conditions are not known. 

• As is intuitively clear, the speed of convergence depends on the distribution of the interac- 
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TABLE III 
Standard deviations for Fig. [8] 



model/pdf 


std for N = 5/10/20/40/100 [dB] 


sumprod/beta 


5.6/3.9/2.7/1.9/1.2 


sumprod/R 


7.6/5.6/4.0/3.0/1.9 


sumprod/L 


6.1/4.2/2.9/2.1/1.3 


prod/beta 


19.7/19.6/19.6/19.6/19.5 


prod/R 


11.7/11.4/11.2/11.0/10.9 


prod/L 


14.2/14.0/13.6/13.8/13.8 



tions as well as its parameters. 
• The fact that the convergence of the quadratic form © to a log-normal random is insensitive 
to the probability law of the individual entries of its component matrices implies that the 
sum-product model is capable of explaining the log-normality of the shadowing phenomenon 
for a wide variety of different channel interactions and radio propagation scenarios. 

VI. Conclusion 

We have studied the applicability of a general sum-product signal model as a basis to explain 
the fluctuation of locally averaged received signal power in mobile communications, i.e., shadow 
fading. It was shown by simulations that as the number of multiplicative "layers of channel 
interactions" (coupling matrices) in the model increases, the distribution of the local average 
power tends to log-normal distribution, as observed in measurements. It was demonstrated that 
only a few layers (e.g. 3 — 5) are needed to produce approximately log-normal power, whereas 
with the conventional model a product of more than 20 random attenuations are needed to produce 
the same goodness-of-fit. Furthermore, the sum-product model does not make the implausible 
assumption that all multipaths experience the same cascade of random attenuations. The sum 
model, which addresses 'local' shadow fading, emerges as a special instance of the general 
model in the case, where the intermediate propagation process does not change and only the 
powers of the plane waves at the receiver end vary. The sum-product model, however, can explain 
log-normality of shadow fading in both local and global scale, and is therefore more general. 
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Fig. 9. Impact of distribution of interaction attenuations to 101og 10 (P) for N = M = 5. The straight dotted line displays the 
normal distribution. Small value of the test statistics indicates good fit. Critical value of the Kolmogorov-Smirnov test for 95% 
significance and 10 3 samples is shown for reference. 



However, more work is needed to derive an analytical justification for the justification of the 
convergence of the sum-product model to a log-normal random variable. 
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